Non-thermal statistics in isolated quantum spin clusters after a series of perturbations 
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We show numerically that a finite isolated cluster of interacting spins 1/2 exhibits a surprising 
non-thermal statistics when subjected to a series of small non-adiabatic perturbations by external 
magnetic field. The resulting occupations of energy eigenstates are significantly higher than the ther- 
mal ones on both the low and the high ends of the energy spectra. This behavior semi-quantitatively 
agrees with the statistics predicted for the so-called "quantum micro-canonical" (QMC) ensemble, 
which includes all possible quantum superpositions with a given energy expectation value. Our 
findings also indicate that the eigenstates of the perturbation operators are generically localized in 
the energy basis of the unperturbed Hamiltonian. This kind of localization possibly protects the 
thermal behavior in the macroscopic limit. 
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The subject of equilibration in completely isolated 
quantum systems has recently seen a significant new 
wave of interest [H-llil, which is motivated by the re- 
alization that the present-day experimental capabili- 
ties are steadily approaching the limit of probing the 
equilibration dynamics of isolated many-body quantum 
svstemsfTBj. In the near future, however, the systems 
available for decoherence-free equilibration experiments 
are likely to be limited to those that have a large num- 
ber of quantum levels but not too large number of par- 
ticles. An isolated system of ten g-bits already fits into 
the above category. Other relevant examples are small 
clusters of nuclear spins, or the systems of a few cold 
atoms. The present paper reveals a surprising statistical 
behavior of such systems under the action of a series of 
small non-adiabatic perturbations. 

If a gas of classical particles in an isolated container is 
subjected to a series of non-adiabatic perturbations, it is 
expected to return to thermal equilibrium once the per- 
turbations stop. The condition for this equilibration is 
that the dynamics of the particles is sufficiently chaotic. 
The quantum counterpart of the above process entails 
one essential complication. Namely, a non-adiabatic 
change of the Hamiltonian of an isolated quantum sys- 
tem generically broadens the participation of the energy 
eigenstates in the wave function of the system. Repeated 
non-adiabatic manipulations can therefore induce a very 
broad participation of the energy eigenstates, which, in 
turn, may lead to a superposition of microcanonical en- 
sembles corresponding to different temperatures [16|. If 
after such manipulations the system is left to itself, the 
occupations of its energy eigenstates and hence the rel- 
ative weights of different microcanonical ensembles will 
not change with time. The small subsystem behavior 
within such a system would be rather unconventional. 

As the above discussion suggests, the appropriate 



choice of the statistical ensemble depends on the prepa- 
ration of a quantum system. For example, energy mea- 
surement, even if inaccurate, is likely to lead to a usual 
microcanonical-like ensemblefsj. However, the ensem- 
bles that might emerge in the course of unitary quan- 
tum evolution associated with non-adiabatic perturba- 
tions have not yet been fully characterized. This line of 
questioning has previously motivated one of us [13, [ill 
to investigate the so-called "quantum micro-canonical" 
(QMC) ensemble [l9l|. which is, in a sense, opposite to 
the conventional microcanonical ensemble. The QMC 
ensemble includes all quantum superpositions in an iso- 
lated quantum system that have a given energy expec- 
tation value i?av, he. the participating eigenstates are 
not limited to the narrow energy window around -Eav, hr 
contrast with the conventional microcanonical ensemble. 
The QMC ensemble is also distinctly different from the 
conventional canonical ensemble. It leads to unconven- 
tional statistics [13, [iM I20l - l26j . which is summarized in 

S3- . 

In this paper, we show that a finite isolated quan- 
tum system with the initial state selected from the 
canonical ensemble permanently departs from the expo- 
nential canonical statistics after a series of small non- 
adiabatic perturbations. The emerging statistics semi- 
quantitatively agrees with the QMC statistics. Our find- 
ings also indicate that the eigenstates of the perturbation 
operators are generically localized in the energy basis of 
the unperturbed Hamiltonian. This peculiar kind of lo- 
calization possibly contributes to protecting the conven- 
tional thermal behavior in the macroscopic limit. Fi- 
nally, we propose that one can obtain the experimental 
evidence of the non-thermal statistics by performing adi- 
abatic magnetization of spin clusters. 

We have performed numerical investigations of the 
periodic chain of 16 interacting spins 1/2 in a time- 
dependent external magnetic field. The Hamiltonian of 
the problem is H = Ho + 'Hp{t), where 
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FIG. 1: Time dependence of the external magnetic field Hx{t). 

and Hpit) = H^{t)ElliSf. Here Sf , Sf , Sf are the 
spin operators for the ith lattice site, the nearest neighbor 
couphng constants are = 2 and Jx = Jy = —1, and 
Hx {t) is a strong time-dependent external magnetic field 
along the a;-direction, which is turned on and off in a 
series of rectangular pulses shown in FiglTJ The above 
choice of the coupling constants is relevant to nuclear 
magnetic resonance experiments in solids (e.g. (28j). 

In general, the energy expectation values of quan- 
tum superpositions increase under non-adiabatic pertur- 
bations. If the perturbations are too strong, then a finite 
quantum system reaches the uniform infinite temperature 
distribution after very few such perturbations with the ef- 
fect of each perturbation being largely non-universal. In 
a search for universal behavior, we have investigated the 
effect of a series of small non-adiabatic perturbations. 

We introduce the non-adiabatic perturbations through 
the time dependence of the magnetic field Hx{t) shown 
in FigHJ The strong magnetic field Hp = 10 acts on the 
system during time interval iom and then it is switched off 
during time toff,n S> t, where t 1 is the characteristic 
timescale of one-spin motion for the Hamiltonian "Ho, and 
n is the pulse number. We used the same ton for each 
pulse but have chosen <off,ri randomly between 100 and 
110. As discussed later, making toff,n constant, would 
make the perturbation sequence periodic, thereby leading 
to the dynamical localization of the occupations of the 
eigenstates under the perturbations. 

In order to illustrate the universal trend towards the 
QMC statistics, we present the results for two qualita- 
tively different kinds of pulse sequences characterized by 
different values of ton- The "short-pulse sequence" has 
very short ton — 0.02 such that Hp ton <C 1, which means 
that during time ton each spin rotates in magnetic field 
Hp by a small angle. In this case, the spread of the occu- 
pations of eigenstates is expected to be of the first order 
in the above angle, while the drift of the average energy 
associated with the Hamiltonian is of the second or- 
der. Therefore, initially, the energy spread is expected to 
grow faster than the energy drift. The second sequence 
contains "mTr-pulses" characterized by Hp ton = w", 
where m is an integer number. In the limit Hp — >■ (X), an 
mTT-pulse simply rotates each spin by rnvr, and hence the 
energy of the Hamiltonian T-Lq remains exactly the same. 
However, since, in reality. Hp is finite, the WTr-pulse in- 
troduces a small perturbation to the system associated 
with the non-commutation between and l-Lp. In com- 
parison, with the short-pulse perturbation, the m7r-pulse 



perturbation is expected to couple more different quan- 
tum states and hence be more effective in forcing the 
system to explore its Hilbert space. The larger m, the 
stronger the perturbation due to the m7r-pulse. Below, 
we use the sequence of 47r-pulses, for which ton — 0.47r. 

We denote the wave function of the system at the end 
of the nth ton-toff cycle as 

k 

where $fc are the eigenfunctions of the system and Cj, 
the complex amplitudes. The corresponding energy ex- 
pectations values e'^^ are 

£;w = 5:£;.pi"\ (3) 

k 

where are the eigenenergies, and p)f' = IC^""*]^ are 
occupation numbers of the respective eigenstates. 

We assume that the spin cluster was initially in con- 
tact with a thermal environment characterized by tem- 
perature Tq ~ 1 {ks = 1) and then was slowly isolated. 
Therefore, for the initial wave function 4'o, we choose 

the complex amplitudes C^."^ — -^e ^^o 6*"*=, where ak 

are random phases, and Z — exp(— i?fe/To). In prin- 
ciple, since the initial thermal state is a mixed one, the 
perturbations-induced ensemble should be obtained by 
averaging over many pure state evolutions with different 
random phases ak- However, given the large number of 
eigenstates, the evolution of a single p ure state already 
contains the required statistics [27j, |29| . The individual 

occupations for a single pure state strongly fluctu- 
ate between adjacent energy levels. In order to obtain 

in) 

the statistics of p\. , we organize the eigenstates into 
energy bins, each containing Nl adjacent energy levels. 
For each bin, we plot p{E), where p = X]fc'"Pfe"' 

^ — ^bin „(re) ■ 

For our pulse sequences, the effect of the nth pertur- 
bation cycle can be expressed as ^„ — C/„5'„_i, where 

Un = exp [-iiUo - Hp J2k Sk)ton] exp(-i'HoWn)- 
Starting with the "thermal" initial wave function ^Pq, we 
calculate the the wave functions ^I^n with the help of the 
complete diagonalization of the Hamiltonians Hq and 
Ho - HpJ2k ^k- 

(n) 

Our numerical results for the statistics of p^. gener- 
ated by the short-pulse sequence are shown in Fig. [21 
Figure [H^a) shows the initial exponential distribution. 
Figure [2Ub) indicates how the deviations from the expo- 
nential statistics begin to develop. Figure [2fc) contains 
a clear case, when the distribution is much better de- 
scribable by the QMC statistics than by the exponential 
statistics. Figure ^d) represents the approach to the in- 
finite temperature state, where all eigenstates are nearly 
equally occupied. 

Figure [3] presents the result of the application of the 
47r-pulse sequence. In comparison with the short-pulse 
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FIG. 2: (Color online) Occupations of the eigenstates of Ho 
after the sequence of n short pulses. Each point represents 
an energy bin of 128 adjacent eigenstates. The QMC and the 
exponential fits are calculated as described in the text. See 
for the semilog version of these plots. 
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FIG. 3: (Color online) Same as Fig. [2] but for the sequence of 
n 47r-pulses. 



sequence, the 47r-pulse sequence generates the QMC dis- 
tribution after much fewer pulses and with more pro- 
nounced deviations from the exponential distribution. 

We have also investigated smaller one-dimensional and 
two-dimensional disordered clusters, where the results 
were qualitatively the same(27|. 

In both Fig. [2] and [Sj we compare the nu- 
merically generated statistics with the the exponen- 
tial statistics p^"\E) = exp{—E/Tn) expected for the 
canonical ensemble and with the QMC statistics 

(")(£;) ^ 1 + Xn{E - Si'i'^) , where both T„ and A„ 
are determined by matching the average energies of the 
respective theoretical ensembles to si^^ — see 

Overall, the numerical results presented in Figs. [2] and 
[3] clearly support the notion that the ensembles gener- 
ated by a series of small non-adiabatic perturbations of 
isolated finite spin systems approach the infinite temper- 
ature distribution through a stage, which is much better 
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describable by the QMC statistics than by the canoni- 
cal exponential statistics. The fact that the 47r-pulse se- 
quence generates less deviations from the QMC statistics 
presumably implies that its effect is closer to the action 
of a random rotation in the Hilbert space. 

Central to understanding the emergence of the QMC 
statistics and the deviations from it is the competition 
between the spread in the occupations of the energy 

(n) 

eigenstates and the drift of the average energy E^v to- 
wards the infinite temperature value. The former reflects 
the degree, to which the quantum superpositions explore 
the "energy shell" in the Hilbert space associated with a 
given value of £:iv' in Eq.©. The QMC ensemble repre- 
sents the limit of random sampling on the above energy 
shell 17|. All possible eigenstates, including those on the 
high energy end of the spectrum, should participate in 
the QMC sampling, but, in order for this to happen, the 
effect of the perturbations should spread sufficiently fast 
towards the high energy eigenstates. 

In a related development presented in Fig. |4j we have 
found that -EI"'' drifts very differently for the periodic 
perturbations with the constant value of the delay time 
toS,n = 100 and for the aperiodic perturbations with the 
random choice of foff.n between 100 and 110. Given that 
the system is expected to relax completely between two 
pulses separated by ioff.n ^ t, both the periodic and the 
aperiodic pulse sequences could have been expected to 
produce the same energy drift. We came to the conclu- 
sion that the observed difference is related to the local- 
ized nature of a typical perturbation operator C/„ in the 
energy eigenbasis of the Hamiltonian "Hq . If the perturba- 
tions are periodic, i.e. Un = Ui, and if the eigenfunctions 
of Ui are localized as suggested, then, for equal to one 
of the eigenstates of Hq, no matter how many perturba- 
tions Ui one applies, the occupations pk would become 
non-zero only within the "energy localization length" for 
the eigenstates of Ui. We believe that the incomplete 
saturation of the energy drift in the periodic case seen in 
Fig. [5] is related to the finite width of the energy spec- 
trum, because, in the localization regime, the transport 
through finite disordered systems is exponentially sup- 
pressed but still possible. For the randomized choice of 
toS,m each operator t/„ is expected to possess a set of 
localized eigenstates, but these eigenstates should be dif- 
ferent for different C/„. Therefore, the effect of a sequence 
of such perturbations involves successive projections be- 
tween the above localized eigenstates which leads to the 
continuous spread of the occupations of the eigenstates 
along the energy axis. The inset of Fig. 0] supports the 
energy localization picture by presenting the eigenstate 

participation ratio / = J2k 1^'^"'!^ • In the case of 
periodic pulses, / nearly stops growing at the same time 
as the energy drift slows down. 

The above localization property is naturally expected 
under the assumption that the effect of each perturba- 
tion Un has a character of disordered one-dimensional 
hopping along the energy axis. The disorder in this case 
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FIG. 4: (Color online) Drift of the average energy i?av for 
the periodic and randomized sequences of short pulses. Inset: 
the participation ratio of the eigenstates (see the text) . 



should originate from the complex structure of the ma- 
trix elements of Un associated with the non-integrability 
of the perturbation problem. A similar kind of dynami- 
cal localizationjsO, Isij has been investigated, e.g., in the 
context of the ionization of Rydberg atoms |33. l33|. 

As the number of spins in the system increases, 
we expect the drift of i?av to scale linearly with Ng , 
while the spread of the eigenstate occupation numbers is 
more likely to have a diffusive character and scale as \/Ws. 



Thus it should be more difficult to abandon the canonical 
ensemble for larger systems. In the large-iV^ limit, the 
eigenstate occupation profiles associated with the above- 
mentioned energy localization may further contribute to 
protecting the conventional thermal behavior j27j. 

If a QMC-like statistics is generated experimentally in 
a cluster of Ns spins 1/2, how can one detect it? A 
possibly effective way of doing it is to perform the adi- 
abatic magnetization of the system with the subsequent 
measurement of the total magnetization. The above pro- 
cedure is implemented numerically in [l^] . 

In conclusion, we have presented numerical evidence 
for emerging QMC-like statistics in isolated spin clusters 
after non-adiabatic perturbations and for the dynamical 
localization property of the perturbation operators. Both 
of the above properties need to be further investigated in 
terms of their dependence on the type and the size of 
quantum systems, on the perturbation routines and also 
in terms of their stability in the presence of decoherence 
and subsequent quantum measurements. 

We are grateful to the bwGRiD computing cluster and 
to the Innovation Fund Frontier at the University of Hei- 
delberg for supporting this work. 



Note added: A very recent preprint of Bunin et al. 34 1 
considers the issue of energy drift vs. eigenstate occupa- 
tions spread in a similar context. 
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SUPPLEMENTARY MATERIAL 

Note: Reference numbers in the text below are from 
the reference hst of the main article. 



A. Quantum micro-canonical ensemble. 

Quantum micro-canonical (QMC) ensemble is defined 
for A^-dimensional Hilbert spaces to include all possible 
superpositions of eigenstates 



N 



(SI) 



fc=l 



with a given energy expectation value i?av, i-e. 



N 
i—k 



(S2) 



"All possible superpositions" means uniform joint proba- 
bility measure on the (2A'^ — 2)-dimensional manifold con- 
strained by condition (IS2p and by the normalization con- 
dition X^ilfc IC'fcP = 1 in the 2A^-dimensional Euclidean 
space of variables {ReCfc,ImCfc}. 

The QMC ensemble allows no exchange of energy or 
particles with the environment. Therefore, at first sight, 
it appears to be similar to the conventional micorcanoni- 
cal ensemble. The crucial difference, however, is that the 
conventional microcanonical ensemble, when applied to 
quantum systems, limits the participation of the energy 
eigenstatates to a narrow energy window of unspecified 
small width (which however, should be large enough to 
include many eigenstates). In contrast, the QMC ensem- 
ble allows all eigenstates to participate as long as the 
low-energy eigenstates balance the high-energy ones to 
satisfy condition ([S2|) . The broad participation of energy 
eigenstates in the QMC ensemble may also appear to be 
similar to the canonical ensemble. Here, however, clear 
quantitative differences exist, which, for example, can be 
seen in Figs. 2 and 3 of the main article. 

The statistical properties of the QMC ensemble is 
more convenient to discuss in terms of the occupation-of- 
eigenstate variables pk = |C/cp and the phases of quan- 
tum amplitudes ak defined by Ck = \Ck\e^"'' ■ The two 
QMC constraints can then be expressed as 



N 



y^^EkPk 



E, 



and 



N 



(S3) 



(S4) 



i=k 



These constraints impose no restriction on phases at, 
which, therefore, can be chosen completely randomly in 



the interval [0, 27r). It can also be shown that the uni- 
form joint probability measure on the QMC manifold in 
the space of variables {ReC/c,ImCfc} translates into the 
uniform joint probability measure in the space of vari- 
ables {pk}- 

Application of the QMC ensemble requires two steps: 
(I) Calculation of the marginal probability distribution 
for the occupation numbers pfc, which we denote as 
Pk {Pk ) , and which gives the corresponding average value 
(pfc); and (II) Calculation of the density matrix for a 
small subsystem. Both steps where done analytically for 
large N in Ref. [17], and the finite- corrections were in- 
troduced in Ref. [18]. The numerical results presented in 
the present work focus on step (I) and indicate that the 
statistics for (pk) generated by a series of non-adiabatic 
perturbations becomes close to QMC. (In the case of 
the conventional canonical and microcanonical ensem- 
bles, the values of (pk) are simply postulated.) Step II is 
discussed at the end of this section. 

The QMC ensemble in the smal\-pk approximation (the 
most general case) leads to the following marginal prob- 
ability distribution[17] 



As a result. 



PkiPk) 



(Pk) - 



-Npk[l+\(Ek-Ee,v)] 



N[l + XiEk - E,^)]- 



(S5) 



(S6) 



Here A is a parameter that can be determined by sub- 
stituting Eq. ((S6|) into the averaged versions of either 
Ea. (IS3l) or Eq. (|S4p and then solving the resulting equa- 
tion for A numerically. It should be noted that the expo- 
nentially decaying distribution PkijPk) has maximum not 
at Pk = {pk) but rather sX pk =0, which means that the 
individual values of pk for any particular superposition 
belonging to the QMC ensemble strongly fluctuate with 
respect to (pk)- For this reason, in analyzing the numer- 
ical results, we use the energy binning procedure, which 
groups together adjacent levels having very close val- 
ues of energies and, therefore, according to Eq. ljSS]) . 
nearly the same probability distributions Pk{pk)- We 
convert formula (jS6l) into the following approximation 
for the total occupation of the quantum states within 
the bin at energy E 



p{E) 



Nl 



N[l + \{E - E^,)] 



(S7) 



and then calculate A numerically by substituting p{E) 
into the the "binned" version of Eq. (jS3p . 

One may be concerned that the integrals / p{E)dE and 
/ Ep{E)dE w\ihp{E) given by Eq. ([S7| appear to diverge 
at large E. In this regard, it is necessary to remember 
that the QMC ensemble is defined for large but finite 
Hilbert spaces which always have the maximum value of 
E terminating the integration range. As the size of the 
system increases, the high energy tail oi p{E) needs to be 
radically balanced by the occupations of the low energy 
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levels to deliver the given value of -Bav It was shown 
in Ref.[17], that, for the systems of macroscopic size, the 
QMC ensemble would lead to a generic condensation into 
the ground state of the system. (See sections II. E and 
II. F of Ref.[17] for more details.) In the present work, our 
numerically generated QMC-like ensemble is far from the 
above condensation regime. 

The perturbations of the external magnetic field in our 
simulations drive the system towards the infinite tem- 
perature equilibrium. In this work, however, we are not 
interested in what happens after infinitely many pertur- 
bations. The question we ask is: what happens if, af- 
ter n perturbations, the system is left completely iso- 
lated. A classical system in this case would thermally 
equilibrate, because it would be describable by a micro- 
canonical ensemble with the value of energy reached after 
the last perturbation. In the quantum case, however, as 
we demonstrated in this work, the resulting ensemble is 
QMC-like for a range of values of n. Since the occupa- 
tions of eigenstates of isolated quantum system do not 
change with time, the resulting QMC-like ensemble will 
not evolve further, which, in turn, implies that the sta- 
tionary density matrices of small subsystems within such 
an isolated system will deviate from the predictions based 
on the conventional thermal equilibrium. 

In this work, we use many small perturbations, which 
change the energy expectation value of a quantum su- 
perposition only by small amount but, at the same time, 
force the system to sample the Hilbert space. We suppose 
that our perturbation sequences generate the QMC-like 
statistics, because the the QMC statistics describes the 
properties of the overwhelmingly typical quantum super- 
position with a given energy expectation value, in very 
much the same way as the Boltzmann-Gibbs distribution 
represents the statistics of an energy shell in the classi- 
cal phase space. The difference between the quantum 
and the classical cases is that in the classical case, the 
phase space trajectory of an isolated system explores the 
energy shell dynamically, while, in the quantum case it 
does not do it by itself - only the perturbations force the 
system to do it, albeit with a simultaneous increase of en- 
ergy. After sufficiently many perturbations and for not 
too large quantum systems, we expect the quantum state 
of our system to become representative of the "Hilbert 
space energy shell" with the appropriate energy expec- 
tation value. A representative sampling of Hilbert space 
energy shells by the technique of small perturbations may 
encounter practical limitations as the number of particles 
in the system increases, because, in this case, the volume 
change of the Hilbert space energy shell after each small 
perturbation should increase exponentially. As indicated 
in the main article, this question requires further inves- 
tigation. 

All characteristics of small subsystems, such as den- 
sity matrix or correlation functions, involve integrals con- 
taining the product v{E)p{E), where v{E) is the den- 
sity of states for the entire system. The Gibbs equilib- 
rium for small subsystems requires that iy{E)p{E) is nar- 



rowly peaked aroimd _Eav[18]. This is always the case for 
the microcanonical ensemble and is also the case for the 
canonical ensemble in the macroscopic limit [18]. At the 
same time, the QMC ensemble in the macroscopic limit 
leads to a double-peak structure of u{E)p{E): one peak 
corresponds to the above-mentioned condensation to the 
ground state, while the other peak corresponds to the 
infinite temperature state associated with the maximum 
of ^{E). (In this work, the energy corresponding to the 
infinite temperature is E = 0.) For the finite systems 
considered in this work, the product of v{E)p{E) is not 
narrowly peaked for the canonical ensemble. Neither it 
exhibits two narrow peaks for the QMC ensemble. In 
such a case, distinguishing the QMC from the canonical 
ensemble at the level of a subsystem is possible but more 
involved because of the non-universal features associated 
with the finite size effects for either of these two ensem- 
bles. We plan to discuss the properties of subsystems 
elsewhere. 

Finally, we would like to remark, that the energy lo- 
calization property reported in this manuscript may lead 
to an exponentially decaying tails of the occupation of 
eigenstates even for the randomized pulse sequences. In 
the macroscopic limit and on the higher-energy end of 
p{E), this exponentially decaying tail would look like 
the canonical ensemble, thereby leading to the narrowly 
peaked shape of v{E)p{E). Such a shape, in turn, would 
justify the Gibbs equilibrium for small subsystems. 

B. Numerical results for small disordered spin 
clusters. 

The Hamiltonians "Ho and 'Hp(t) considered in the 
main article are translationally invariant. In addition, 
in the quantization basis of operators S^, neither %q nor 
'Hp{t) has oH'-diagonal elements coupling the quantum 
states having even and odd numbers of spins with pos- 
itive a;-projections. As a result, the quantum basis can 
be decomposed into 32 non-mixing subspaces, each asso- 
ciated with one of the 16 possible values of the inverse 
lattice wave vectors and with the even/odd property. Be- 
low we break both of the above symmetries by introduc- 
ing additional disorder. In comparison with the 16-spin 
chain considered in the main article, the disordered sys- 
tems investigated below have smaller numbers of spins 
but larger sizes of mixing Hilbert spaces. These exam- 
ples are included also to illustrate that the Bethe-ansatz 
intergrability of the Hamiltonian T-Lq for spin chains is not 
a determining factor in the emergence of the QMC-like 
statistics. 

In this section, we modify the Hamiltonian 'Ho as fol- 
lows: 

NN 

i<j i 

In comparison with Eq.(l) of the main article, the lat- 
tice indices in the interaction part are modified to ac- 
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FIG. SI: Occupations of the eigenstates of Ho for a 13-spin 
chain with disorder after a sequence of n short pulses. Each 
point represents an energy bin of 16 adjacent eigenstates with 
additional averaging over 48 independent choices of the initial 
random phases Ok- The QMC and the exponential fits are 
calculated as described in the main article. 



commodate the nearest neighbor interaction (NN) for an 
arbitrary lattice dimensionality, but the main difference 
is the last term. It contains random local field along 
the z-direction. In the context of NMR experiments, this 
term may represent the disorder in the values of chemical 
shifts. All other parameters of the simulations remain the 
same as in the main article unless indicated. [Note: In 
the presence of the above disorder, only m7r-pulses with 
even m perfectly conserve the average energy associated 
with T-lo in the limit Hp oo.] 

In Figs. lSll and lS2[ we present the results for the short- 
pulse and the 47r-pulse sequences, respectively, acting on 
a periodic 13-spin with local fields hf chosen randomly 
in the interval [—0.2,0.2]. In the both figures, the re- 
sults are additionally averaged over 48 independent ran- 
dom choices of the phases for the initial wave func- 
tion 5*0. This additional averaging reduces the statistical 
noise, but as mentioned in the main article, the energy- 
binned result for a single randomly chosen set of {a^} 
already exhibits the same average occupations of eigen- 
states. 

In Fig. IS3| we present the results for the sequence of 
27r-pulses acting on a periodically closed 4x3 square 
lattice with local fields chosen randomly in the interval 
[-0.1,0.1]. 



C. Series of 47r pulses for periodic spin chains of 
different length 

In Fig. 3 of the main article, we present the statistics 
for the occupations of eigenstates generated by different 
numbers of 47r-pulses on the 16-spin periodic spin chain 
with anisotropic interaction. Here we fix the number of 
47r-pulses (n = 5) while changing the number of spins 
Ns in the spin chain from 10 to 16. The initial ensemble 
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FIG. S2: Same as Fig. ISll but for a sequence of n 47r-pulses. 
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FIG. S3: Occupations of the eigenstates of Ho for a 4 x 3 
square lattice with disorder after a sequence of n 27r-pulses. 
Each point represents an energy bin of 32 adjacent eigen- 
states. The QMC and the exponential fits are calculated as 
described in the main article. 



in each case is chosen in the same way as in the main 
article — it corresponds to the thermal occupations of 
the eigenstates with temperature T = I. The results are 
presented in Fig. IS41 



D. Adiabatic magnetization 

Here we describe how one can detect the QMC statis- 
tics using the adiabatic magnetization routine mentioned 
in the main article. 

If the magnetic field is increased sufficiently slowly, 
the eigenstates evolve, their energy change, but their oc- 
cupations do not change much. At a large enough value 
of the magnetic field, the energy of each eigenstate is 
mainly due to the Zeeman energy of the total magnetiza- 
tion in the magnetic field. By adiabatically magnetizing 
a single cluster, and then doing the measurement of its 
total polarization = J2i ^ one can obtain Fm^ , the 



8 








M 

X 




-10 -5 5 10 -10 -5 5 1 



E E 

FIG. S4: Occupations of the eigenstates of Ho for periodic 
spin chains of different length A'^a after a sequence of five in- 
pulses. Each blue point represents one energy bin of adjacent 
eigenstates. The eigenstates in each case were divided into 
512 equal energy bins. The QMC and the exponential fits are 
calculated as described in the main article. 



probability for the occurrence of different values of M^. 
The number of the quantum states with a given value 
is = (i/2jv,-Af^)T°(i/2jv,+M^)! ' Therefore, the oc- 
cupation of each of these states is Pm^ = Fm^ / ■ The 
values of Pm^ are expected to exhibit nearly the same 
statistics as the occupations of the cluster eigenstates 
before to the magnetization begins. 

We have performed a numerical imitation of the above 
procedure. We started from the numerically generated 
QMC-like state of the 13-spin chain shown in Fig. 182^0). 
The magnetization was performed by increasing the mag- 
netic field Hx from to 10 in steps of 0.1 separated by 
the waiting time 1000. The result is shown in Fig. [S5] 
The signature of the QMC statistics in such an im- 
perfect adiabatic procedure is rather clear. One ex- 
pects Pm^ — exp(Ma;_ffa;/r'), if the magnetization pro- 
cess starts from the canonical distribution, and Pm^ ^ 
[1 — X'Hx{Mx — Mx,av)]~^, if the magnetization starts 
from the QMC ensemble. Here M^r.av = J2m ^m^M^ 
evaluated for the numerically generated ensemble. For 
the respective fits shown in Fig. IS5[ both T' and A' were 
determined numerically from the simulations-generated 

value of Mr^av 

The inset of Fig. [S5] shows the probability Fm^ of mea- 
suring a given value of ■ This plot indicates that the 
measurement outcomes would be dominated by the states 
around = due to the large degeneracy dM^ of these 
states. Therefore, a large number of single cluster mea- 
surements should be made in order to extract the sig- 
natures of the QMC statistics around the minimum and 
the maximum values of M^- If many equivalent elec- 
tronic spin clusters are available for simultaneous ma- 
nipulations, then one may try to bring nuclear spins in 
contact with these clusters and then measure the distri- 
bution of Mx through its effect on the NMR lineshape. 



FIG. S5: Occupations of individual quantum states (blue cir- 
cles) with a given value of spin cluster polarization af- 
ter numerical imitation of adiabatic magnetization for a 13- 
spin chain with small disorder. The adiabatic magnetiza- 
tion process started from the eigenstate occupations shown 
in Fig. IS2f c'). The green triangles and the red squares repre- 
sent, respectively, the exponential and the QMC fits obtained 
as described in the text. The inset shows the probability of 
measuring a given value of for the simulation results (blue 
dots), exponential fit (green triangles) and the QMC fit (red 
squares) . 



E. Semi-logarithmic plots corresponding to Figs. 2 
and 3 of the main article and to Figs. [STl [S2] and [S3] 
of the supplement. 
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FIG. S6: Semilog version of the plots shown in Fig. 2 of the 
main article. The plots represent the results of the short-pulse 
sequence acting on a 16-spin chain without disorder. 
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FIG. S7: Semilog version of the plots shown in Fig. 3 of the 
main article. The plots represent the results of the 47r-pulse 
sequence acting on a 16-spin chain without disorder. 



FIG. SIO: Semilog version of the plots shown in Fig. IS3I The 
plots represent the results of the 27r-pulse sequence acting on 
a 4 X 3 square spin lattice with disorder. 




FIG. S8: Semilog version of the plots shown in Fig. ISII The 
plots represent the results of the short-pulse sequence acting 
on a 13-spin chain with disorder. 




FIG. S9: Semilog version of the plots shown in Fig. IS2I The 
plots represent the results of the 47r-pulse sequence acting on 
a 13-spin chain with disorder. 



